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Abstract 

The equation of state of symmetric nuclear matter is addressed starting both 
from a realistic interaction derived from nucleon-nucleon scattering processes 
and from a low-momentum effective potential. The approach is based on finite 
^sO , temperature Green's functions. The internal energy per particle is estimated 

from the summation of diagrams and through the Galitskii-Koltun's sum rule. 

<N 
-t— > . 

■ At present, there is not a unique and reliable theory of nuclear interactions. 

From the study of nucleon-nucleon scattering processes different models for the 
basic interaction between nucleons have been developed, but difficulties arise 
, when they are used in calculations which aim to reproduce the properties of 

bound nuclei. These models all reproduce the experimental N-N phase shifts 
up to the threshold for pion production with high accuracy, but differ in the 
short- wavelength region, where less experimental constraints are available. 
, In deriving the properties of dense nuclear matter, if we start from a bare 

N-N potential we have to take into account the complex correlations induced 
by strong interactions. In parallel with the development of diagram summation 
techniques, a variety of effective potentials have been proposed. Apart from 
f--^ | the traditional phenomenological N-N forces, to be used in a mean-field ap- 

proximation, such as the Skyrme and the Gogny interactions, in the last years 
renormalized low- momentum potentials have been introduced pQ. They are re- 
stricted to a subspace of the Hilbert space, the so-called model space, in which 
they incorporate the short-range correlations between nucleons. If the cutoff is 
sufficiently low, they turn to be independent of the starting N-N interaction. 

In this work we present a self-consistent scheme based on finite temperature 
Green's functions and apply it to study the equation of state of symmetric nu- 
clear matter. Wc first consider a bare nucleon-nucleon interaction, namely the 
; ^ ■ CD-Bonn potential, and compare the results with the case in which a renormal- 

ized potential Viow-k HI is used. 

It is possible to construct a consistent approximation starting from a suitably 
chosen generating functional 3 . Such approximation schemes automatically ful- 
fill thermodynamic relations, including the Hugenholz-Van Hove and Luttinger 
identities. For nuclear interactions the generating functional must at least in- 
clude ladder- type diagrams: this choice leads to the in- medium T- matrix, the 
approximation scheme adopted in our study. 

The in-medium two-particle scattering matrix T is defined as (for simplicity 
we skip spin and isospin indices, as well as energy and momentum dependence): 

T = V + VG% C T . (1) 
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Here V is the interaction potential and GVf is the non correlated two-particle 
Green's function, just constructed as a product of two dressed single-particle 
propagators. 

The T-matrix accounts for the multiple scattering between nucleons and 
thus permits to resum the short range correlations induced in the dense system. 
We remark that the use of the dressed propagator implies the presence of a 
nontrivial dispersive self-energy which leads to a broad spectral function (5J- 
The full two-particle propagator is then approximated as follows 

g 2 = G% c + Gl c T G% c + exchange terms . (2) 

Eqs. (JIJ and J2J, together with the Dyson equation, form a scheme in which all 
ingredients are calculated iteratively: the scattering matrix, the single-particle 
self-energy, which is expressed in terms of the scattering matrix, and the single- 
particle propagator. 

Most of the properties of the system can be derived from these quantities. 
The (total) internal energy per particle can be calculated as the expectation 
value of the Hamiltonian H = Hkm + H pot . The two terms read (we now make 
explicit the momentum and energy dependences) 



m \-v f d3p dcJ p2 



A(p,w)f(w), (3) 



where V represents the volume, A(p,u>) is the spectral function and f{io) the 
Fermi-Dirac distribution, and 

V f d 3 P d 3 k dn 
J (2tt)3 (2,)3^ b(0) 

xIm{((k|T*(P,fi)|k) - (k|T fl (P,fi)|-k)) G? fl (P,k,ft)} , (4) 

where Q and P are the total energy and momentum of the pair of nucleons, k 
the exchanged momenta and 6(0) the Bose-Einstein distribution. 

An alternative and simpler way to determine the energy is the Galitskii- 
Koltun's sum rule 00 



E 1 f d 3 p did [ p 2 

' ■ A(p,u)f(u) , (5) 
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which for conserving approximations is equivalent to the direct calculation 
©+©■ 111 the presence of three-body forces, however, this is not valid and 
the expectation value of the Hamiltonian has to be computed to estimate the 
internal energy of the correlated system. 

The ladder expansion can be as well derived from a generating functional 
<!>[G, V], from which pressure and entropy at finite temperature can be obtained 

We performed calculations for the internal energy per particle in two ways: 
from the Galitskii-Koltun's sum rule (J5J and from diagram summation, i.e. eqs. 
© and (0J. For the two methods we employ both a realistic nucleon-nucleon 
interaction, the CD-Bonn potential, and the effective interaction Vi ow -k from 
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Figure 1: Diagrammatic expression for the expectation value of the potential 
energy in the case of the full T-matrix approach (upper figure, the sum runs 
over the number of interaction lines n in the diagram) and for a second order 
approximation (lower figure). 



|2] . In the first case we consider the whole expansion of the propagator in ladder- 
type diagrams, which lead to the diagrammatic expression for the interaction 
energy that appears in Fig. H(up). The energy depends on dressed propagators 
(solid lines in the figure) and on the two-body potential (wavy lines). In the 
case of the effective interaction we only include diagrams up to the second order 
(Born approximation), which are shown in Fig. ^(down). Since Vi ow -k already 
accounts for the correlations in the dense system the higher order diagrams do 
not yield significant contributions and can be neglected [2J. 

We restrict ourselves to symmetric nuclear matter at zero temperature, 
though we remark that this formalism can be applied to finite temperature 
systems as well. Results for the internal energy per particle as a function of 
density are shown in Fig. in which calculations from the diagram summation 
are compared to the Galitskii-Koltun's sum rule. The calculations for the two 
potentials are in agreement up to 0.8po- At higher densities, the low-momentum 
calculations give too attractive energies and differ substantially from the full T- 
matrix ones, exhibiting no minimum. On the contrary, the full calculation with 
a realistic nucleon-nucleon potential shows a saturation point, which however 
does not coincide with the experimental value (as expected since three-body 
forces are not included.) 

Concluding, we have studied the equation of state of symmetric nuclear 
matter both with the two-body CD-Bonn potential and with a renormalized 
interaction Viow-k- When the effective potential is used, the energy behavior 
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Figure 2: The internal energy per particle at zero temperature as a function 
of the density (in units of the empirical saturation density po = 0.16 fm -3 ). 
The solid lines represent the expectation value of the Hamiltonian the 
dashed lines are the results obtained from the sum rule lp|l. 



shows that Vi ow ~k alone at high densities does not correctly account for all 
correlations. This would not happen if a larger cutoff were chosen, even if in 
this case Viow-k would depend on the starting nucleon-nucleon realistic force. 
This dependence, which can be amplified or suppressed by changing the cutoff, 
can possibly help in testing the reliability of the different models of the bare 
nucleon-nucleon interaction. 
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